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ABSTRACT 

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear 
to correlate with the physical properties of their hosts, including the mass of the 
dark-matter halos. At higher redshifts, we observe the growth of SMBHs indirectly 
through the identification of high redshift quasars. However information on their hosts 
is more difficult to obtain. In this paper we determine the masses of the halos that 
host the high redshift quasars (at z > 4) by comparing the rate of growth of quasar 
density with that predicted by the Press-Schechter mass function. The host mass 
determined depends on how the ratio between SMBH and host halo mass evolves 
with redshift. Under the assumption that the ratio between SMBH and halo mass 
does not evolve with redshift, we find a host halo mass of M — iQii ''^'' ^. Even if the 
quasars shine at their Eddington limit, this host mass is significantly smaller than that 
seen at lower redshifts in the local universe. Indeed we find that the null-hypothesis, 
of a constant ratio between SMBH and halo mass at all redshifts, can be ruled out 
at greater than a 5-sigma level. SMBHs must therefore have contributed a larger 
fraction to the host mass in the past. This finding is consistent with expectations from 
models of self limiting SMBH growth. When we include the redshift evolution of the 
ratio between SMBH and halo mass, we find larger halo masses of M ~ j^q12.4±o.3^ 
in combination with a ratio between SMBH and host halo mass that increases with 
redshift in proportion to ~ (1 -I- z)^-^ are required to be consistent with both local and 
high redshift observations. We also investigate the restrictions placed on the critical 
linear overdensity of quasar hosts at their epoch of virialisation and find that it cannot 
exceed the traditional value of Sc = 1.69 by more than a factor of two. Finally, we find 
that the high redshift quasars are hosted by fluctuations on scales that have a variance 
of (SM/AI) =2 — 3, corresponding to (3 — 4.5)-sigma fluctuations in the density field. 
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1 INTRODUCTION 

Tlie Sloan Digital Sky Survey has discovered luminous 
quasars at redshifts as high as 2 ~ 6.4, i.e., when the 
universe was only a billion years old. The super-massive 
black-holes (SMBH) powering these quasars have been es- 
timated to have a mass of about 10^ solar masses. How- 
ever questions regarding the galaxies that host these high 
redshift quasars have remained largely unanswered. To un- 
derstand the formation and evolution of quasars and the 
super-massive black-holes that power them, one needs to 
determine several important physical parameters (such as 
the quasar lifetime, the ratio of black-hole mass to halo 
mass, and the efficiency and rate of accretion during the 



luminous phase), as well as the evolution of these parame- 
ters with time. Attempts to answer these questions gener- 
ally consider the quasar luminosity function, as it provides a 
tracer of the density of quasars with different luminosities as 
function of cosmic epoch (e.g. Haehnelt, Natarajan & Rees 
1998; Haiman & Loeb 1998; Kauffmann & Haehnelt 2000; 
Volonteri, Haardt & Madau 2003; Wyithe & Loeb 2003). 
However, all these analyses are model dependent and im- 
plicitly assume either (a) the quasar lifetime and its evo- 
lution with redshift, and/or (b) the form, normalisation 
and evolution of a relation between SMBH mass and the 
characteristic velocity of the host galaxy. The exceptions 
are studies at low redshift that utilise the quasar two-point 
correlation function (Martini & Weinberg 2001; Haiman & 
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Hui 2001; Croom et al. 2004). At high redshift dynamical es- 
timates have been made in a few individual cases (Barkana 
& Loeb 2003; Bertoldi ot al. 2003). 

Locally, direct estimates of SMBH and host mass 
can be made through observations of galaxy dynamics. 
These observations reveal a correlation between SMBH 
mass and the characteristic velocity of the surrounding stel- 
lar spheroid (e.g. Mcrritt & Fcrrarcsc 2001; Tromaino ct 
al. 2002), and by extension of the host dark matter halo 
(Ferrarese 2003). These characteristic velocities determine 
the dynamical mass, so that there axe also corresponding 
correlations between SMBH and host mass. Any proposed 
scenario for SMBH evohition must reproduce this behaviour 
and hence the correlations provide important clues regard- 
ing the physics of formation of SMBHs. 

However at higher redshift, objects collapse out of a 
denser back-ground. The characteristic velocity of a viri- 
alised object of a given halo mass is therefore larger if it 
formed at higher redshift. This begs the question of whether 
the fundamental correlation is between SMBH mass and the 
hosts characteristic velocity, or between SMBH mass and 
dynamical host mass. Obviously knowledge of which cor- 
relation is fundamental is critical for our understanding of 
the astrophysics of SMBH evolution. Unfortunately since 
these dynamical observations can only be made for relatively 
nearby galaxies, this question cannot be resolved via direct 
observation. 

Here, within the paradigm of standard concordance cos- 
mology, we show that by associating halo mass to quasar lu- 
minosity, the formation rate of luminous quasars in the high 
redshift universe can be used to constrain the mass of the 
dark matter halos that host them. Wo thus provide a frame- 
work within which one can attempt to answer the question 
of whether host mass or velocity is the determining factor in 
the evolution of a SMBH. We believe this procedure holds 
significant promise for the future when the observations im- 
prove. 

In § 2 we begin with the null-hypothesis that the ratio 
between the masses of SMBH and the halo does not evolve 
with redshift. Under this assumption, we find the host mass 
to be about 10"-^*°-^Mo. Estimates for SMBH mass pow- 
ering the luminous SDSS quasars yield ~ W^Mq assuming 
output at the Eddington limit (the Eddington limit pro- 
vides a lower limit on SMBH mass, and hence a lower limit 
on the mass ratio). The resulting SMBH to halo mass ratio 
is therefore much larger for the SDSS quasars than for lo- 
cal galaxies. Indeed, the null hypothesis can be rejected at 
a significance greater 5-sigma, i.e. SMBHs in the past con- 
tributed a larger fraction of galaxy mass than those today. 
This is one main conclusion of the paper. Next, we allow for 
the SMBH to halo mass ratio to evolve with redshift (§ 2.6). 
In this more general case we estimate a halo mass for high 
redshift quasar hosts, M ^ 2q12.4±o.3^ g^^j that the ra- 
tio between SMBH and host halo mass should increase with 
redshift as ~ (1 -|- z)^'^ in order to be consistent with ex- 
trapolation from local observations. Estimates for the high 
redshift quasar lifetime are discussed in § 3. We also inves- 
tigate variance of the linear power-spectrum on the scale 
of density fiuctuations corresponding to the high redshift 
quasar hosts in § 4. Finally, we discuss the the implications 
of the rate of high redshift quasar density evolution for the 



value of the critical linear overdensity at host virialisation 
in § 5. Some concluding discussion is given in § 6. 

Throughout the paper we adopt the set of cosmolog- 
ical parameters determined by the Wilkinson Microwave 
Anisotropy Probe (WMAP, Spergcl ct al. 2003), namely mass 
density parameters of Qm = 0.27 in matter, O;, = 0.044 in 
baryons, Qa = 0.73 in a cosmological constant, and a Hub- 
ble constant of Ho = 71 kms~^ Mpc~^. For the primordial 
power-spectrum of density fluctuations, we adopt a power- 
law slope n — 1, and the fitting formula to the exact transfer 
function of Bardeen et al. (1986). It turns out that our re- 
sults are most sensitive to the cosmological parameter as, 
which is the amplitude of the linearly extrapolated power- 
spectrum on scales of Sh^^Mpc. We present results for dif- 
ferent erg wherever appropriate to illustrate the range of this 
dependence. 



2 EVOLUTION IN QUASAR DENSITY 

The Press-Schechter (1974) mass function (with the mod- 
ification of Sheth & Tormen (2002) that will be adopted 

throughout our discussion) yields the number density A'^(> 
M{z),z) of dark matter halos above some mass M{z) at red- 
shift z. If luminous quasars reside in a fraction e of such dark- 
matter halos, then the observed number density of quasars 
is given by the product of two factors; N{> AI{z),z) and 
r = e mm{tq/ H^'^ (z) , 1} where tq is the (unknown) quasar 
lifetime and H~^{z) is the Hubble time (see e.g., Estathiou 
& Rees 1998). As a measure of the rate at which luminous 
quasars appear, we use the logarithmic derivative {B) of 
tN{> M{z), z), defined as 

dlogr d\ogN{> M,z) 7 d In iV(>M,g) 
~ dz (1-Fa) dlnM ^ ' 

Here we have assumed that the halo mass varies as M oc 
(l-^-zy at fixed luminosity M1450. This choice is convenient 
for our discussion because observations do reveal an expo- 
nential decline in the quasar population with redshift sug- 
gesting approximate constancy of B. Our analysis, of course, 
relies on the applicability of the Press-Schechtor mass func- 
tion (as modified by Sheth and Torman 2002). It hcis been 
shown that this analytic formalism provides an excellent de- 
scription of the halo mass function found from numerical 
simulations (Jenkins et al. 2001). In particular, if the nu- 
merical mass function is expressed in the appropriate vari- 
ables then it is independent of epoch, which is a defining 
feature of the Press-Schechter formalism. Evaluation of B 
using the Press-Schechter mass function therefore provides 
an accurate description of the redshift evolution of massive 
dark-matter halos in a ACDM cosmology. 

The mass function A'^(> M, z) leads to a steeper slope B 
as M is increased. This is illustrated by the grey dot-dashed 
lines in the left panel of figure 1 which show the evolution 
of A''(> M,z) (arbitrarily normalised at 2 = 4.8) for masses 
of 1O^°M0, 10" Mq, IO^^Mq and 10^^ Mq. These curves 
should be compared to the observed evolution in the den- 
sity of luminous (M1450 < —26.7) quasars between redshifts 
of 2 ~ 3.7 and 2; ~ 6 from the Sloan Digital Sky Survey 
(Fan et al. 2001;2003;2004), which is also summarised in the 
same panel. Here we are assuming that the identification of 
luminosity M14B0 with halo mass M does not vary (that is. 
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Figure 1. Constraints on the evolution of the density of high-rcdshift quasars. Left: The density of quasars with M1450 < —26.7 as 
a function of redshift (Fan et al. 2001;2003;2004). The solid line shows the best fit exponential deelinc *(< 26.7, z) oc 10^°baXz 
the full data set, with Bobs = —0.49. For eomparison the evolution in the density of halos (normalised to unity at z = 4.8) with 
masses of IO^'^Mq, IO^^Mq, IO^^Mq and 10^'^ Mq are shown by the grey dot-dashed curves in order of decreasing steepness. Right: The 
corresponding a-posteriori differential (grey lines; left axis) and cumulative (black lines; right axis) probability distributions for Bobs- 



7 = 0). This forms the null hypothesis in this work. The ef- 
fects of relaxing this assumption are discussed later in § 2.6. 

The curves corresponding to 10^^ Mq and 10^^ Mq have 
logarithmic slopes that lie at the extremes of the range al- 
lowed by the data. It is therefore clear that if we can mea- 
sure the exponential slope in the formation of high redshift 
quasars, then we can determine the mass of their host dark- 
matter halos. By using only the logarithmic slope B, we have 
removed the dependence on the absolute value of quasar life- 
time. However, B does depend on the form of the redshift 
evolution of r (and on 7). While e could also change with 
z due to various effects (dust obscuration, beaming angle 
etc.), we expect the dominant additional contribution to B 
to come from tq. The z dependence of tq can be handled by 
using two physically motivated forms for the evolution — 
which bracket the reasonable range of possibilities and 
can be parameterised by t cx (1 -|- z)" with ^ a ^ 3/2. 
First, if the quasar lifetime is determined by the mass e- 
fold timescalc of the SMBHs, then tq is independent of red- 
shift, T oc 1/H^^{z) and a ~ 3/2. Second, if the quasar 
lifetime is determined by the dynamical timescale at z, then 
tq Ri H~^{z) making t independent of redshift and a ~ 0. 
(This is also true if tq > H~''{z).) 



2.1 Evolution constraints 

The evolution shown in figure 1 is well fitted by an expo- 
nential decline (Fan ct al. 2001) of the form 

*(Mi450 < -26.7, z) cx 10^°'>=^". (2) 

The right hand panel of figure 1 shows the a-posteriori dif- 
ferential (grey curves) and cumulative (black curves) proba- 
bility distributions for the observed exponential slope Sobs- 
These distributions were computed as follows. For each value 
of Bobs, we find the normalisation which maximises the 
product of probabilities from each redshift bin (assuming 
Gaussian error bars). This product represents the likelihood 
for Bobs. A flat prior probability for Bobs was then assumed. 



allowing calculation of a-posteriori distributions for Bobs- 
Distributions were estimated using (i) the whole data set 
and two subsets of the data; (ii) data with z < 5, and 
(Hi) data with z > 4.5. We find that Bobs ~ -0.49 ± 0.07, 
Bobs ~ -0.52 ± 0.15 and Bobs ~ -0.53 ± 0.20 describe the 
evolution within the full data set and two subsets respec- 
tively showing internal consistency. The best fit is also plot- 
ted in the top left panel to guide the eye. 



2.2 Mass constraints 

From equation (1) there is a one-to-one monotonic corre- 
spondence between halo mass M and exponential slope B. 
The a-postcriori probability distributions for host halo mass 
A4 may therefore be found by noting that 

dP dB dP 

dM°^ dM B=B„, -"^ ' ^ ' 



where (dB/dM) was computed using equation (1) at the 

central redshift within the luminosity function data. The dif- 
ferential (grey lines) and cumulative (dark lines) a-posteriori 
distributions for M are plotted in the left hand panel of fig- 
ure 2. If = 3/2 (solid lines), the observed evolution in 
the density of bright quasars implies host halos with masses 
of W^^-^^°'^Mq. If a = (dashed lines), then the halo 
mass estimates and their upper bounds are about a fac- 
tor of 2 smaller; we get lO"-®*"-^M0. Thus under the null- 
hypothesis of a fixed 7 = 0, the exponential slope of the 
high redshift quasar luminosity function leads to the deter- 
mination of masses of halos that host high redshift quasars 
to within a factor of a few. In particular, we stress that the 
result does not rely on any a-priori assumptions about the 
relation between quasar luminosity and SMBH mass, about 
the relation between SMBH and halo mass, or about the 
halo density profile. 

As an aside, we note that an alternative approach 
to equations (3) [and (5)] for calculation of dP/dM [and 
dP/d{etq)] would be to choose prior probabilities for M and 
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Figure 2. Constraints on the mass of tlie halo that hosts a high redshift quasar. Left: The a-posteriori differential (grey lines) and 
cumulative (dark lines) probability distributions for M (solid lines a = 3/2, dashed lines a = 0; 7 = 0). The upper axis shows the 
corresponding values for fraction of halo mass contributed by a lO^AfQ black-hole. Right: The corresponding differential (grey lines) and 
cumulative (dark lines) probability distributions for etq. All curves in this figure were evaluated for its = 0.84. 



etq that are fiat in the logarithm, and likelihoods for M and 
etq based on the comparison of the corresponding slope with 
the data. We find that this approach gives nearly identical 
results to those presented in this paper. 



2.3 SMBH to halo mass ratio 

We now explore several further consequences. For the esti- 
mation of SMBH mass, it is usually assumed that quasar 
emission is isotropic, and that emission is at the Edding- 
ton rate, resulting in SMBHs powering the highest red- 
shift quasars having inferred masses (Fan et al. 2001) of 
~ 10® Mq. These estimates are consistent with dynamical 
estimates based on emission line profiles (Willott, McLure 
& Jarvis 2003). We have therefore labeled the upper axis 
of the left panel in figure 2 with the fraction of halo mass 
contributed by a 10® M© black-hole, allowing the curves in 
this panel to represent the a-posteriori probability distri- 
butions for this fraction as well as M. We find that the 
full data set implies SMBHs contribute a fraction of about 
10"^' 



^(Mbh/lO^Mo) and ~ 10" 



HMbh/10®Mo) of 



the halo mass for a = 3/2 and a = respectively. These 
fractions are larger than those found by Ferrarese (2003) for 
local ~ IQ^'^Mq galaxies, which are M^h/M = lO""' '^, lO"'^ 
and 10~^'^ respectively under the assumptions of singular 
isothermal halos, Navarro, Frenk & White (1997) (NFW) 
halos, and halo masses derived from galaxy-galaxy lensing 
(Seljak 2002). 

This result has two possible interpretations. First, the 
SMBHs are accreting at well above their Eddington rate 
and/or with high efficiency so that the black-hole mass has 
been overestimated (this disagrees with dynamical studies, 
Willott et al. 2003). A second, more agreeable, interpreta- 
tion is that SMBHs may contribute a larger fraction of the 
halo mass at higher redshifts. Indeed this latter scenario 
is independently supported by observations of quasar host 
galaxies at z ~ 2 which suggest black-hole masses that are 



significantly larger with respect to their hosts than in the 
case of local galax:ies (Groom et al. 2004). 



2.4 Sensitivity to as 

Of the observable cosmological parameters, the relationship 
between B and M is most sensitive to as. To illustrate the 
extent of this dependence we have repeated our analysis us- 
ing values of as = 0.76 and (78=0.92. (These values bound 

the 2-CT range of the best fit constraints derived from WMAP 
plus large-scale structure and Lya forest data Spergel et 
al 2003.) We find that the constraints on the quasar host 
halo mass vary by a factor of ~ 4 within the 2 — a range for 
Us (see figure 3). 



2.5 Marginalised distributions 

One can also marginalise over systematic uncertainty in 

dP/dM due to a and as; 



dP 
dM 



<.3/2 <.o 

/ da I 

Jo Jo 



das exp 



-(erg - 0.84)2 



2(0.04)2 



dP 
dM 



(a, as), (4) 



where {dP/dM){a,as) was determined using equation (3), 
and we have used a fiat prior probability for a in the range 
< a < 3/2 combined with a Gaussian probability distri- 
bution for CTg. We find (left hand panel of figure 4) a value 
of M = 10" '^'''°-^Mo, which is our best estimate for the 
halo mass of the high redshift galaxies that host the quasars 
under the null hypothesis of 7 = 0. Assuming a 10® Mq cen- 
tral black-hole, corresponding to accretion at the Eddington 
limit (Willott, McLure & Jarvis 2003), this corresponds to a 
black-hole to halo mass ratio of lO"^ '^*^'^, which is inconsis- 
tent with local estimates (Ferrarese 2003; ~ 10~® assuming 
an NFW profile) at greater than 7 — cr. If luminous quasars 
shine near their Eddington rate over a range of redshifts, we 
therefore conclude that the SMBH to halo mass ratio must 
increase with redshift. 
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Figure 3. Constraints on tlie mass of dark matter lialos wliicli liost quasars. The a-posteriori diflferential (grey lines) and cumulative 
(dark lines) probability distributions for M obtained using values of as =0.76 and as = 0.92. These correspond to the 2 — cr range for 
IT8 determined from WMAP (Spergel ct al. 2003). The solid and dashed curves in these panels correspond to a = 3/2 and a = 0. In each 
panel the upper axis shows the corresponding values for the fraction of halo mass contributed by a 10® Mq black- hole. 
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Figure 4. Constraints on the mass of dark matter hales which host quasars. Left: A-posteriori probability distributions for M that 
have been marginalised over a and as (see equation 4). The upper axis shows the corresponding values for the fraction of halo mass 
contributed by a 10® Mq black-hole. Right: The 2.5, 16, 50, 84 and 97.5 percentiles of the cumulative marginalised probability for mass 
M as a function of 7. The grey regions represent masses outside the l-cr range for each 7. For comparison, we show curves (dashed lines) 
representing quasar host halo mass at high redshift derived assuming the local black-hole halo mass ratio (Ferrarese 2003; for a lO^-^M© 
and 10® M0 SMBH, and the NFW profile) plus an evolution of this ratio with redshift that is proportional to (1 -|- z) >. 



2.6 Mass estimates assuming an evolving SMBH 
to halo mass ratio 

In evaluating the derivative B, we have thus far assumed a 
luminosity M1450 to be associated with a fixed halo mass 
(7 = 0). However if the halo mass M housing quasars of lu- 
minosity M1450 varies with redshift, then the evaluation of 
B must include the extra term in equation (1) of the form 
7M(2:)/(1 + z) X dlogio A''(> M, z)/dM (where we have as- 
sumed that the halo mass varies as M cx (1 + z)'' at fixed 
luminosity M1450.) In the right hand panel of figure 4 we 
show the 2.5, 16, 50, 84 and 97.5 percentiles of the cumu- 
lative marginalised distribution (equation 4) for halo mass 



as a function of 7. The grey regions represent the masses 
outside the l-cr range for each 7. Smaller values of 7 lead 

to larger estimates of the mass (Note that the y-axis is re- 
versed with 7 decreasing from bottom to top). If the evolu- 
tion of halo mass M housing a fixed black-hole mass follows 
M oc {1 + zy between the local and high redshift universe 
then we can estimate M at high redshift as a function of 7, 
by extrapolating the local relation (Ferrarese 2003) between 
halo and black-hole mass. The resulting curves (taking the 
case of an NFW profile) are plotted for black-hole masses 
of lO* ''''Af0 and 10® Mq (dashed Unes from left to right). By 
comparing the two constraints we find 7 ^ —(1.5 — 2), lead- 
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Figure 5. Constraints on the mass of the halo that hosts a high redshift quasar. Left: The a-posteriori differential (grey lines) and 
cumulative (dark lines) probability distributions for M (solid lines a = 3/2, dashed lines a = 0; 7 = —1.5). The upper axis shows the 
corresponding values for fraction of halo mass contributed by a lO^AfQ black-hole. Right: The corresponding differential (grey lines) and 
cumulative (dark lines) probability distributions for etq. All curves in this figure were evaluated for its = 0.84. 



ing to estimates of halo mass that are ~ 3 — 6 times larger 
than our evolution-free (7 = 0) estimate. Note that these 
larger masses do not weaken our result that SMBHs com- 
prised a larger fraction of galaxy mass at high redshift, as 
this behaviour is explicit when 7 < 0. Interestingly, values of 
7 < follow naturally from models where SMBH growth is 
self limiting through feedback on galactic gas (e.g. Haehnelt, 
Natarajan & Rees 1998; Wyithe & Loeb 2003). 

We see from Figure 4 that if the evolution of the ra- 
tio between SMBH and halo mass (or more correctly, the 
ratio between quasar luminosity and halo mass) can be de- 
scribed as a powerlaw in redshift, then an extrapolation from 
local observations combined with our evolution analysis im- 
plies a value of 7 ~ —3/2. In the left panel of Figure 5 
we plot the differential and cumulative probability distri- 
butions for M assuming 7 = —3/2. As mentioned in the 
previous paragraph we find larger masses than were derived 
in Figure 2 under the assumption of non-evolving ratio; we 



get M = 
(a = 0). 



10 



12.B±0.2 



Mq (a = 3/2) and M = 10 



12.3±0.2 



Mr. 



3 QUASAR LIFETIME 

If a fraction e of dark matter halos contain SMBHs, then 
the total lifetime of the quasar can be estimated (Martini & 

Weinberg 2001; Haiman & Hui 2001) by dividing the quasar 
number density >J'(Afi45o < —26.7, z) by e times the rmmber 
density A'^(> M) of halos larger than M, and then multi- 
plying by the Hubble time (for < H~^). The a-posteriori 
probability for the product etq is 



dP 
d(etci) 



dM J d^ ' 



dB 
dM 



dP 



-B— Sobs 



,(5) 



where (dP/d'i/) is the observed Gaussian probability for 
>]/. Under the assumption of a SMBH to halo mass ratio 
that does not evolve with redshift (see Figure 2, right hand 
panel), we find lifetimes of lo^-sio-^e-i and lo^-SiO-^e-i 



years (respectively) for a = 3/2 and a = 0, in a cosmol- 
ogy where as = 0.84. In addition (not shown) we have com- 
puted distributions for eiq corresponding to ag = 0.76 and 
(78=0.92. We find constraints that vary by a factor of 2 rel- 
ative to the case of as = 0.84. These results imply that 
if all dark matter halos contained SMBHs at high redshift 
(e = 1) then the preferred quasar lifetime of (10* - 10^) 
years is significantly shorter than both the Salpeter time 
[about 4 X 10'^(eoff /O.l)/;^"*" years] for accretion at the max- 
imal rate (?? = 1) with tes = 10% efficiency of conversion 
from mass to energy, as well as estimates of the quasar life- 
time at lower redshifts [(lO" - 10*) years, see Martini (2003) 
for a summary] . The small value of etq might therefore indi- 
cate that not all dark-matter halos at high redshift contain 
supermassive black-holes, in contrast to the situation locally 
(Kormendy & Richstone 1995). Indeed, if the high redshift 
quasar lifetime were (10^ — 10*) years, as seems to be the 
case (Yu & Tremaine 2002) at 2 ~ 2, this would imply that 
only 1 in (10 — 10"^) dark matter halos at 2; > 4 contained 
a SMBH. Alternatively, there could be a larger number of 
obscured quasars at high redshift, giving the impression of 
a smaller quasar lifetime. 

On the other hand, we have suggested that a situation 
where the SMBH to halo mass ratio does not evolve with 
redshift will be inconsistent with local observations. As a 
result it should not be surprising that the lifetime derived 
when 7 = is inconsistent with other observations. In § 2.6 
we derived results for the halo mass under the assumption 
that the SMBH to halo mass ratio evolve eis a power-law 
in redshift [oc (1 -I- z) '] and found that larger masses of 
M ~ 10^^'''='=°-^Mq are obtained for 7 = -1.5. We have also 
derived the a-posteriori probability distribution for quasar 
lifetime in this case (right panel of Figure 5), and found that 
the larger halo mass leads to a longer inferred lifetime. We 
find tq = io5-9±0'4e-i years (a = 3/2) and tq = lo^'SiO-S^-i 
years (a = 0). For occupation fractions of unity (e = 1) these 
lifetimes are marginally consistent with, though still smaller 
than the lifetime of 10^ — 10* years inferred at lower redshift. 
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4 THE VARIANCE OF DENSITY 

FLUCTUATIONS CORRESPONDING TO 
HIGH REDSHIFT QUASAR HOSTS 

An alternative way of studying the hosts of the high red- 
shift quasars is to compare the evolution of their density 
with the evolution of the density of halos that correspond 
to scales with a fixed variance in the lineaxly extrapo- 
lated power-spectrum. In Figure 6 we show curves of den- 
sity [tN{> M{a),z)] against redshift for values of variance 
a = 2, 2.5 and 3. Since we are interested in the devia- 
tion of the observed results from the theoretical curves for a 
constant a, we have normalised both the theoretical curves 
(from Press-Schochtcr) and the data to unity at an inter- 
mediate redshift z — 4.8, and thus express the evolution in 
terms of a dimensionless density parameter [f2(z)/f2(4.8)]. 
The left and right panels correspond to values of a = 3/2 
and a = respectively. We find that the values of a = 2 — 3 
bracket the range of evolution, implying that for a linearly 
extrapolated critical overdcnsity of Sc ^ 1.69(1-1-2), the high 
redshift quasar dark-matter host galaxies formed from ~ 3- 
sigma to ~ 4.5-sigma density fluctuations. The rareness of 
these fluctuations is consistent with the large masses (well 
in excess of the non-linear mass-scale) which were inferred 
in § 2.6. 



5 THE LINEAR OVERDENSITY FOR 
QUASAR HOSTS 

When calculating the Press-Schechter mass function, we 
have so far adopted the conventional value for the linear 
overdcnsity at halo virialisation of 5c = 1.69, which is ap- 
propriate for a spherical collapse at the time when shells 
evolve to zero radius. This value of 5c, which corresponds 
to a non-linear overdensity of about 178, agrees with re- 
sults from N-body simulations (Jenkins et al 2001). How- 
ever, in general, non-spherical top-hat models virialisc at 
a time when the linear overdensity reaches a value greater 
(Engineer, Kanekar & Padmanabhan 2000) than 5c = 1.69. 
As we showed in § 4, luminous quasars are rare systems 



forming from greater than 3-sigma fluctuations. It is there- 
fore possible that they form only in hosts of unusual overden- 
sity. In figure 7 we have plotted contours of i? as a function 
of 5c and M. [for the case of a = 3/2 and as = 0.84, and £is- 
suming 7 = (left panel) and 7 = —3/2 (right panel).] The 
shaded grey regions show the excluded values of 5c < 1.69, 
and B (68% range). The most likely value of the variance 
is plotted as a thick dashed line. We see that 5c cannot be 
too difi^erent from 1.69, or else the host dark matter halos 
would be unacceptably small. This indicates that collapse 
to an unusually overdonsc halo is not taking place in halos 
that host SMBH formation and quasar activity. An absolute 
upper limit on 5c can be obtained by noting that if the mass 
of the SMBH is 10® M©, then host masses must be larger 
than about 10^° M©. In this case, values of Jc ^ 3 are not 
allowed. If the SMBH wore restricted to contain less than 
10% of the gas component, then 5c ^ 2. 



6 DISCUSSION 

In this paper we have estimated the mass of high redshift 
quasar host dark matter halos by comparing the rate of 
quasar evolution with the Prcss-Sclicchtcr mass function. We 
have found that in the case of the null hypothesis where the 
SMBH to halo mass ratio does not change with redshift, that 
the implied halo mass is M = 1O" ''='=°-^M0. This mass is 
significantly smaller than the mass of local halos that house 
a 10^ Mq SMBH, the mass befieved to be powering the SDSS 
quasars. Indeed our results rule out the null-hypothesis at 
a significance greater than 5— a. We therefore conclude that 
the SMBH to halo mass ratio must increase towards higher 
redshift, i.e. SMBHs contained a larger fraction of the host 
galaxy mass at earlier times. 

Having demonstrated that SMBHs at high redshift must 
have contained a larger fraction of the host mass than 
SMBHs observed today, wc allowed the SMBH to halo 
mass ratio to vary with redshift. In this case it is possi- 
ble to achieve consistency between observations of high red- 
shift quasars and local SMBH masses. We find that these 
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Figure 7. Contours of the exponential slope _B as a function of the linear overdensity for quasars 5c, q and host halo mass M (we assumed 
a = 3/2). The shaded regions are excluded by the requirement that 5c, q > 1.69, and by the limits on -Bobs derived from the redshift 
evolution in the quasar luminosity function. The dashed line shows the most likely value for the slope, and the dotted line corresponds 
to 5c,q = 1.69. Left: curves for 7 = 0. Right: curves for 7 = —1.5. All curves in this figure were evaluated for as = 0.84. 



combined constraints imply high redshift quasar host halo 
masses of M = ioi2.4±o.3^ ^j^j^ ^ SMBH to halo mass ratio 
that varies with redshift approximately as (1 + z)^^"^. 

A scenario where SMBHs formed at high redshift con- 
tain a greater fraction of the host galaxies mass is consistent 
with models of SMBH evolution in which SMBH growth is 
limited by feedback during the quasar phase (e.g. Silk & 
Rees 1998; Haehnelt, Natarajan & Rees 1998; Wyithe & 
Loeb 2003) . These models predict a relation between SMBH 
mass and the characteristic velocity of the host which is red- 
shift independent. As a result feedback regulated growth of 
SMBHs leads naturally to a ratio of SMBH to halo mass that 
increases with redshift. Our results therefore support feed- 
back regulated schemes where SMBH growth is dominated 
by accretion during the luminous quasar phase. 

Upon completion, the Sloan Digital Sky Survey will 
have identified much larger numbers of high redshift quasars 
than are currently published . The more accurate luminosity 
functions which will be available should then allow a similar 
and more precise analysis using the model-independent tech- 
nique introduced in this work. Such an analysis may allow 
the variation of host mass with redshift to be determined 
directly for the high redshift quasars. 
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